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This paper is a partial study of currents induced on circular, conducting cylinders by 
narrow radiating slots. First, a brief and general formulation of the radiation fields of 
slots on cylinders is made. Then, the problem of an infinite axial slot is examined thoroughly 
for all cylinder sizes. An expansion for the fields, very close to the slot, on large radius 
cylinders, is also obtained. Sample computations are made, for various ranges of cylinder 
radius, and the order of the errors is discussed. 

The problem of a circumferential slot, of constant excitation, is also considered. An 
asymptotic expansion obtained for this case yields the surface current distribution for values 
of axial distances that are smaller than the square of the circumference of the cylinder. 

Since one of the objectives of this study is to determine mutual coupling between two 
slots on a cylinder, the last section presents a formulation of the equivalent network in 
terms of the surface and feed line currents. 

1, Introduction 

The problem of determining the interaction of sources with each other or with other 
discontinuities, on curved conducting surfaces, is basically a problem of the fields of the sources 
on the conducting surface. In a previous report Held and Hasserjian [1958] discussed the 
results of an experimental study showing the effect of the tip of a cone with a slot. The con- 
clusions of that study led us to consider the problem of the fields of a slot on a circular cylinder. 

There were several reasons that led us to this approach. Our study of the cone indicated 
that the cone problem was basically a study of fields on curved conducting surfaces with the 
cone tip as a singularity. In adchtion, the mathematical analysis of surface fields on the cone 
was tliought to be a formidable problem. Therefore, a simpler geometrical model was sought 
which could be studied both theoretically and experimentally and then the results extended 
to not only the cone but other geometries as well. Therefore, it is believed that a complete 
study of a circular cylinder will provide fundamental answers to various pi'oblems of curved 
conducting surfaces. 

The formal solution of the fields due to sources on cylindrical structures has been treated 
in the literature [Silver and Saunders, 1950; Sensiper, 1953 and 1957; Wait, 1959]. However, 
expressions for numerical computations have been limited to the determination of the fields 
at large distances from the cylinder surface. Several authors have treated the problem of cur- 
rents on the surface of cylinders [Papas and King, 1949; Lucke, 1953; Nishida, 1960] for infinite 
axial slots. However, their results do not indicate the magnitude of the errors involved in 
the approximations that they employ. 

In this report we provide expressions, for the surface currents due to infinite axial slots, 
for all ranges of cylinder radius. The expressions that are used are all as^miptotic which pro- 
vide estimates for the errors. Comparisons are made between the harmonic series and residue 
series evaluations. 
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The problem of the currents excited by a narrow circumferential slot of constant excitation, 
running entirely around the cylinder, is also treated. The results are again asymptotic expan- 
sions and in this case are limited to axial distances small compared to the circumference of 
the cylinder. 

2. Formulation of the Cylinder Problem 

In this section we shall formulate the expressions for the fields due to slots located on 
conducting circular cylinders of infinite length. Such a formulation has been performed by 
various authors [Silver and Saunders, 1950; Wait, 1959] and is outlined briefly here. 

The general solution, for the field components, may be derived from two scalar functions 
which individually satisfy the wave equation. We shall identify these functions by ^](r,<^,Z) 
and \l/2{T,ip^Z). The coordinate system and the cylinder orientation is shown in figure 2.1. 
The formal solution for the two functions may be written 

^TT J —03 n= — oo 

where an appropriate path should be chosen for the Fourier inversion. The coefficients of 
the series, A^i and A^^j may be determined if the field distribution on a closed surface is known. 
In this problem, the tangential components of the electric field on the surface of the conducting 
cylinder will be zero everywhere except across the slot. If the field across the slot is known, 
then the field distribution on the entire surface of the cylinder is defined, which makes it pos- 
sible to determine the unknown coefficients of the general solution. 

Considering two slot orientations, longitudinal and transverse, and assuming the tangeJi- 
tial electric fields to have only components transverse to the slot, we have the following bound- 
ary conditions on the surface of the cylinder (r—a). For a longitudinal slot 

£'^=0, for all Zand <p 

E^=M<p)f2(Z), for Z,<Z<Z2, <P,<<P<cp2 (2) 

E^=Oj everywhere else. 
For a transverse slot 

Ez=gi(^)g2(Z), for Z(<Z<Z2, ip[<(p<(P2 

Ez=0, everywhere else (3) 

E^=0, for all Z and (p. 

These provide sufficient information to determine the coefficients of the general solution and 
the field components derivable from them. 

For a longitudinal slot, the coefficients of the series of 1/^2 are zero because Ez=0. The 
field components as a function of \p] are : 



^^=i(S^+^0^-^-o 



H^ = — ^ ^ry ) E^= ^— (4) 

cojLte rO(poZ e Or 

'" ajjue brdZ ' er dip 



Matching (1) and (2) by the relation E^= ^7 we find that 
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Figure 2.1. Coordinates and orientation of the 
cylinder. 




^n 



5„,(r) 



where 



^T^ JZi J<Pi 

Therefore, tor the longitudinal slot 



(5) 
(6) 

(7) 



Once (7) is evaluated, we can then evaluate the other field components from the rela- 
tionships in (4). 

For a transverse slot, both scalar functions i/'i and \l/2 have non-zero coefiicients. 
The field components as a function of ipi and 1/^2 are : 

i d^ipu 1 ^^2 



Th 



(jo/jLe rdifbZ iJL dr 
i d''\l/u , 1 ^^2 



TT __ 

^ (jOjjLe drdZ ' jjL rdcp 



(N) 






e dr icofie rdipbZ 

1 bh, , 1 C)V2 



e rd(p ioofjie drbZ 
Using eq (8) and the boundary conditions (3) we obtain from the Ez component 



-^»2 " 



(20)^6)^^2(0 



wher( 



Also, from E^ we have 






A„l T 



«i //f(aVF-nA2 
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(9) 



(10) 



(11) 



or 



n^BAO 



_/€\ HJ^„2(.n 



f) 



Therefore, for the scalar function ypi t of the transverse slot, we have 
and 



(F-n^'Wf (dV^'-n 



f/f 






(12) 



(13) 



(14) 



Substitution of (13) and (14) in (8) will jield the desired field components at any point 
around the cylinder. 

In the expressions of (7), (13), and (14) we have integrals of the form 



h{Z) = J_°'^ e*f^/„(f)Hf (r VF- f)dt 



(15) 



where the path of integration should be specified around the singularities of the integrand. 
These singularities are two branch points located at ^=^k. 

For the Hankel function of the second kind, it is necessary to require that over the path of 
integration of (15) 

-(TT+eXargVF^^e. (16) 

Furthermore, in order to have outward radial propagation and to keep the fields bounded at 
infinity, it is necessary to require that for k complex 



Re [V^'-n>o, im [V^'-n<o. 



For k real, (17) reduces to 



argV^'-f'=0, for IfKA: 
argVF^^=-|for|fi>Z:. 



(17) 



(18) 



Therefore, in discussions of the integral in (15), we shall denote the path of integration by 
C], as shown in figure 2.2. The branch cuts and the path of integration conform with the con- 
straints specified by (16), (17), and (18). Wlien k is real, the path of integration should avoid 
the singularity as shown by the dotted curve in the figure. 

This completes the formulation of the fields of slots on conducting cylinders. In the 
sections that follow, we shall evaluate the magnetic field on the surface of the cylinder for a 
uniformly excited longitudinal slot of infinite length and a circumferential slot running entirely 
around the cylinder. 



C=-k 



C-PLANE 



Figure 2.2. Path of integration in the ^-plane. 



rc=k 
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3. Infinite Axial Slot With Constant Excitation 

A longitudinal slot, of narrow width, can be said to have a constant electric field across the 
slot. Therefore, on the surface of the cylinder 

E,=^^[U{<p+h)-U{<p-fo)]j°'jiZ')S{Z-Z')dZ' (1) 

where Vq is the voltage across the slot at some reference point, 2\l/o is the angular width, a is the 
radius of the cylinder, /(Z) is the axial variation of the voltage along the slot and the functions 
inside the brackets are the usual Heaviside, step functions. Since the integral implies the 
summation of infinitesimal slots, for the purpose of formulating the problem, we shall first 
assume 

/(Z)=5(Z-Z'). (2) 

Therefore, considering (1) and (2) and evaluating the coefficients Bn\{^) defined in (2—6), 
we have 

Substituting (3) in (2-7), we have 

In order to evaluate the fields due to any slot with an excitation /(Z) along its axis, we 
need to integrate the product oi j{Z') and (4) with respect to Z\ 

To obtain the appropriate expression for the longitudinal slot of infinite extent, of 2ai/'o 
width and of constant excitation along the Z-axis, we integrate eq (4) with respect to Z^ from 
— c» to +00. Since 

^J"_^e'f'^-^'>(/Z' = -5(f) (5) 

then (4) becomes 

^n(r,^)-+^^j^^2^^ ^^^ e ^<2).^,^)- (6) 

From this result expressions for any of the field components can be obtained. However, 
we shall concern ourselves only with the numerical evaluation of the magnetic field on the 
surface of the cylinder. 

In order to obtain expressions which are amenable to computations, it is necessary to con- 
sider approximations valid in different ranges of cjdinder radius. We will obtain approxima- 
tions valid in the ranges ka <C<Cl, ka<25 and ka^ 10. 

3.1. Small Radius Cylinders 

The scalar function (6) evaluated on the surface of the cylinder is 

{ N _ — Vq^ r -^0^^ (kg) ^ sin nxl/p cos ncp Hj,^^ (ka) 1 
yp, ^a,^j -— ^ \jkaH['' {ka) hi~^^ ka H'i''{ka)\ 

for n> 1 and Z:a<<l 
Therefore, 
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(7) 



"f^)— (^)[i+ot(^«)^"- («) 



The sum of the infinite series in (9) 



Sr{<p) = 12 



sin nx/zQ cos Uip 



is a symmetric function of (p and is evaluated in appendix I. Therefore, 

using (2 to 7) and (11), the magnetic field on the surface of the cylinder becomes 



H(ff 






-sm'J 



(10) 



(11) 



(12) 



where for \1/o<(p<t 



S,(^) = 1 +(^^HL^ _l^ In (2 si.i(^/2) -v^o cos(^/2)) 
-(^^^+1) 1^(2 sin(W2) +^0 cos(^/2)) . 



(13) 



Appendix I also includes the evaluation of Si for — iAo< v<^o, making it possible to evaluate the 
magnetic field across the slot. Figure 3.1 A shows the plot of sample computations for the 
magnetic field distribution on the surface of cylinders of ka^<il and i^o=l/40, and figure 
3.1B is an expanded plot of the magnetic field in the region of the slot for various values of 
slot widths 24'o- Both figures represent the plot of the bracketed term in (12). 



4''ii<p)-- 



m'ijca) 



~2kaH['\ka) 



-SM. 
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Figure 3.1A. Plot Q/^i= 2kaHi(2>(ka) "^^^^'^^' 



-\.b% -1.0% -0^% 0.5% \0% \.b% 2.0% 

<t> - ANGULAR DISTANCE FROM CENTER OF SLOT 

Figure 3. IB. The expanded plot of \p'x{<p) around slot 
/or ka=0.1. 



3.2. Intermediate Range Cylinders 

For cylinders having radii 0</:a<25 (this also includes the first range), the numerical 
evaluation is obtained by approximating the Hankel functions for large orders and using the 
approximation for the higher indexed terms of the series. A paper by C. H. Papas and R. King 
[1949] shows this evaluation without indicating the order of the errors involved. In this 
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section we attempt to bridge this gap. Furthermore, the results obtained in this analysis 
can be employed to determine the magnetic field across the slot of widths up to 30° with great 
ease, making it possil)le to evaluate the slot impedance. The examples, considered in the 
analysis, are assumed to have a constant electric field across the slot. However, the analysis 
can be applied to other field distributions, provided the axial distribution is constant. 
Taking the expression in (7) and letting ka=p, we write 

For n\ l—(-j >>1 and -=cosh «, we have (see [Watson, 1952] page 243) 

Hr(.)^-i^'-^ .^'«- ^-.. »' [l+O (j^)] (15) 

Therefore, the ratio of the two I'ujictions above becomes 

This leads to an approximate expression for the series in (14), which can now be expressed 



as 



^^ sin 7i\l/Q cos 7i(p Hi;\p) ^ sin nipo cos rup p^ ^ sin n\l/Q cos Uip ^, 



where, for sufficiently large M 



(17) 



It'.<±'^+,^±'~^- (18) 

n=M 'f^ n = M ''^ 

Considering, now, the series 

_ ^ sin n\p() cos Uip 

appearing in (17), we write 

1 ^ sin^cos^ ,, . 

where, for sufficiently large M 

In order to evaluate the upper bounds of (18) and (20), we employ the Euler-Maclaurin 
sum formula. The members of the series in (18) and (20) have the form 

Un)='^- (21) 

Theret'ore, with Bp and (pp{t) as the Bernoulli constants and polynomials respectively, we have 
from Whittaker and Watson [1952, p. 128], 
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where 



00 /»°° 1 m-1 ( ^\V'R 

T..iM=\ JMdx+y,{U)+l^ ^ '>J^^ frKM)+R^ (22) 

n=M JM ^ p = l \^V)' 

Considering (21) in (22), it is immediately apparent that the dominant term of (22) will be 

r r°° cos ifX , 

Jm X 

Integrating this by parts, we have 

J. _ sin M(p.k r^'sinJi^ , 

Furthermore, 

— TTTT- dx< -,^^ , for MiD small 






^riT^f^^^v^la^g^'- 



Therefore, for a sufficiently large M, we can write 

where by the sym.bol, < >, we mean the bounds described above. 
This leads to the conclusion that, for a sufficiently large M 

±;-^=0[<^>} (23) 

We can now^ use (23) to estimate the upper bounds of (18) and (20). Therefore, the 
expression in (17) becomes 



«-I '^^ cos „, [ JPgj+lH-^, im-b'+m *]+«» <^*' 



where 



<-, ^ sm nxpQ cos n(p 
^1=2—1 ; 



(^ _^ sm nxj/Q cos Uip 



The infinite series Si and 5^2 can be expressed in a closed form and are evaluated in appendix I. -< 

The dominant term in the expression for Rm depends on the values of M, <p, and p. 
Therefore, the magnetic field on the surface of the cylinder is 
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where 



M-l, 






r m>{lca) 1 J_/ka\n 



foi" \po<ip<Tr, and 



S2 = ^-^^ + V+-A^^ 



6 2^ 



12' 



The order of the error of (25) can be estmiated from (24). Figures 3.2, 3.3, 3.4, and 3.5 
are sample computations of the bracketed term in (25), which are identified by xpi. These 
computations were performed by a digital computer and the following conditions were 
imposed. 

M=10forO<^a<3 
and 

M=3[ka] ioY3<ka<25. 

By [ka] we mean the integer following ka. 

The results in this section show^ clearly tlie orders of the a])|)roxiinatio]is us(m1 in computing 
the magnetic field on the surface of the cylindei" by using finite tenns of tlw Jiai-monic series. 
Even though one exampk^ of the slot angular width (21/^0=1/20) was considered, the expressions 
are valid for widths up to \I/q=1/2. In addition, for values of ka large, fairly accurate values 
for the surface fields near and across the slot can be obtained. 

These results can also be compared with the residue series evaluations sinc(» computations 
are not restricted to small values of ka. Such a comparison is ina(h' in section 3.4. 

In computing the examples in the figures, the Hankel functions for S^ (see eq (25)) were 
computed by using the recursion fonnulas. For the rcjd part of the Haid^el function, Jn{p), 
the iterative steps were started at index M and M—1, and for the imaginary part of the 
Hankel function, F,j(p),the iteration was begun from n=0, 1. This approach i'(m1uccs possilde 
cumulative erroj's due to the iteration process. 
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Figure 3.3. Real part of the magnetic field on the 
Figure 3.2. Imaginary part of the magnetic field on cylinder surface, for ka = 0.5, 1.0, S.O, 10, \I/q= 1/4-0. 

the cylinder surface, for ka = 0.5, 1.0, 5.0, 10, 

^0 = 1/40. 
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Figure 3.4. 
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cylinder surface. 
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Figure 3.5 A. Phase of the magnetic field (ka = O.S, 
1.0, 5.0). 
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Figure 3.5B. Phase of the magnetic field (ka = loj 



3.3. Large Radius Cylinders 

For cylinders having radii in the range ka> 10, the standard Watson Transformation 
[Bremmer, 1949; Sensiper, 1957; Wait, 1959] can be used quite adequately. The expression 



ypi{a,ip)-- 



Foe y. sinTMAo ^.,^^ HS\ka) 



2'Kka 



n\f^o 



Hrika) 



is a residue series of the contour integral 



hia,<p)-- 



Fpe ri r sin #0 g-^'"''-'" Hj'^ika) 
2wka\_2i J(y' v^'o sin »7r Hl^^'{ka) 



dv\ 



(26) 



(27) 



where C^ is specified in the z?-plane as in figure 3.6. The series in (26) represents the sum of the 
residues of the poles of the integrand in (27) that occur on the real axis. However, the contour 
C[ can be deformed into a new contour Cg, which encloses the zeroes of Hf^ (ka) providing 
another residue series. Figure 3.6 also shows the deformed contour C2 which encloses the 
poles, Vm. The poles v^ are the set of roots of 

m''\ka)=0 
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Figure 3.6. The contours of C[ mid C2 in the v-plane. 





identified by the index m=l, 2, . . . . This new series converges rapidly for ^ and ka large. 
However, it will be necessary to evaluate the integral (27), by another method, for small 
values of <p. 

We first take the residue series about the poles Vm- Since 

and if the zeroes of HJ^^ (ka) are v^ in the lower half plane, then the poles in the upper half 
plane are Vjn=—Vm. Therefore, the integral in (27) expressed as a residue series is 



^'^«''')=-£ 






+ S 



sin Vm^o e 



-iv^(<p—T) 



H^'\ka) 



'^m^O 



Sin... (|-Z/f'(te) 



F( 






sin Vm^Po cos Vm{(p—7r) 



ka ^i Vm^Q sin Vrnir 



m''{ka) ■ 



bv 



(28) 



Note the sign change due to the clockwise rotation of the contour C2. Equation (28) expresses 
the fields on the cylinder surface as attenuated standing waves in the ^-direction. It is also 
possible to express the fields in terms of traveling waves. We take 



sm vw ^ 



(29) 



as a function of the complex variable v. The first sign represents the value on the upper 
half plane and the second sign represents the value on the lower half plane. Considering 
the two parts of (28) and identifying ^f and ^^ as the two sums representing the residues in 
the upper and lower half planes respectively, we have 



ka p=om^l ^m^O 



^+ = 



iV^ 



— ^ S S J?l ^^^ ^^^^ e-''rn [2P^+<P] 

ka p=o w=i Vm^o 



(30a) 
(30b) 



where i?i represents the ratio containing the Hankel functions. 
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The Vm^, which are roots of Hl^^\ka)=0, fall in the fourth quadrant of the «;-plane. 
Therefore, one could say that (30b) represents the waves traveling in the positive cp direction 
and that the integers, p, represent the waves of mode m that have encircled the cylinder p- 
times. For large cylinders the ^=0 term would be sufficient to represent the field since waves 
that have completely gone around the cylinder will be highly attenuated. 

To evaluate i?i, the ratio containing the Hankel functions, we use the approximation for 
the Hankel function discussed in appendix II. 

We need to evaluate 



Ri- 



where the vjs are the roots of Hi;P'(ka)=0. 

From (11) and (12) of appendix II we have 



(31) 



and 



where 

.=z+(|)"'<. 

If we let Z=ka and substitute the above approximations into (31), we have 

-.-(i)"""^ii =-(f)-^ 

at v=v^ 

where the t^s are the roots of co^(t)=0 corresponding to the vjs by the above relationships. 
From the differential equation of the function cx){t), discussed in appendix II, we have 

c^''{t)=to:{t). (33) 

Substitution of (33) in (32) leads to 

/y\2i-i 1 

Substituting this result in (30), »« for —v^ and taking tiie p=0 term we have 



^t=+i(^)(^4)"'±'^' 



e~ 



where 

q^2/3 iir 



(35b) 



and x„ represents the mth root of 
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It is apparent that the expression in (35a) represents the waves travehng- in the negative ip 
direction and in (35h), those travehng in the positive <p direction. 
Therefore, for hirge cyhnders, the total magnetic field is 

Foj- sufficiently large values of (p the serums is rapidly convergent. For instance, with k(i = 25, 
the m = 4 term is less than 1 percent of the m = 3 term for (^=20°. 

To evaluate the integral in (27) for cp small, we again consider the traveling waves to be 
highly attenuated at the end of a 360 degree path. Therefore, we can again substitute 

2ie-'''' for 



sin vir 



in eq (27) and consider the integral along the lower half of the contour. This gives the expres- 
sion foi' the traveling wave along the positive cp direction. The expression now becomes 

We evaluate the ratio of the Haidvcl functions above by using again the first term of the 
asymptotic expansion discussed in appendix 11. We have from (11-11) and (11-12) 

2?rM.-Jp(|)'^(o 



^""'^^^^^i^G^F'^'^^^ 



which gives for v=ka-\-l -^) t 



Hrika) fkaV'^coifl 



This approximation is based on the assumption that / remains finite as (ka) increases. The 
integral in (37) imphes an integration for 0<|/Koo and the approximation in (38) is valid for 
finite values of /. This apparent contradiction is compatible due to the exponentially decaying 
nature of the integrand in (37) and Watson's Lemma. Therc^fore, the major coiitribution to the 
integral is due to the finite values of t. 'I'hus, taking the approximation of (38) and the first 
term of the expansion 

sin Vm^o sin kaxpo f, 1/2 V^^ , , . i , \ 

+~ (^)'" t^ (l-karPo cot fo^„+Mo)!y . . .J (39) 

The integral in (37) with the new variable of integration, t, becomes 

V,e sin ka<I.U kaY\aa. r'"^\-'Ciy" ^' ^!^ ^t (40) 

^'~ 2^ka kaio \2j ' J-w-^ ^'(t) ^ ^ 

We now make the following change of variables : 

x=ka(p=2Tr times the distance in wavelength from the center of the slot 
2xo—2ka\l/o=^2T times the slot width in wavelength 
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c=2 ( — 1 ^^^ = a curvature term 
5 =—i^=: change of variable of integration. 
The expression in (40) now becomes 



where 



,+^^Vo sin Xq eJ2 j 

r fx 1 f-^'+^- £, co(^5) , 
\c 27^^ J-.^-tco w'(^«) 



(41) 



(42) 



The singularities of the integrand above lie on the left hand side of the line of integration 
in the s-plane. Since the origin is not a singularity we can express (42) as a Bromwhich integral 
with T=x/c 



Therefore, with (/' representing the Laplace operator 






We next attempt to find an expansion of 



^ 00 {is) 
in inverse powers of s to obtain an expansion of /(r) for r small. We have 



(43) 



(44) 



(45) 



'it)=^le-'h-t)m%[_l{-ty"J 



Therefore 



ZT(2) 

ns)=4^— 

^^—IS TT{2) 
J^2/S 






I i-isy 



(46) 



(47) 



In order to develop the Hankel functions in their asymptotic expansions, we need to assure 
that the phase of the argument lies in the range — 27r to w. We note that in (43) s will range 
in value from +2oo to —i^. Therefore, we use the identity 



in (47) and obtain 



lfi?{f3[-V^|s''^] 



Fis) = 



Vi^|?^[-V^|.s'^'^] 



(48) 



(49) 



Using the asymptotic expansions of the Hankel function we obtain for their ratio 
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-''"\}+lz-ik^-'ik^--]- (^«) 



- 2 

If we now let Z=-'^Ji - s^^^ and substitute (50) in (44) we obtain 



(51) 



/w=.-.fvv^-.[j=4:+3«,+^.+...]. 

Then employing the basic relation 
it follows that 



where 

X kacp 



(53) 



r= 



c ^/toy 



<■ 



2/ 



Substituting (53) in (41) and letting y= — j-^ — j we finally have 



^,-+i2ilsi^.^.-„.-.,.[(,_|,+^^+. ..)_,(|„_Jj,.+^„+ . . .)] 



2 Xq 



(54) 



The curvature term, the series appearing in (54), is rapidly convergent for small y. It 
approaches unity as (ka) approaches infinity. Therefore, the expression 

which appears in (54) as a multiplying factor to the curvature term, should represent the 
field due to the slot on a flat sheet. However, takmg the same limit of the original expression, 
eq (37), yields the following results: 

where Xc=ka\l/Q and x^=ka(p\ 

We have from Sensiper [1957] or Nishida [1960], with x=ka(p 



Therefore (56) becomes 



lim -. e-'^'^--'' i,Ai^a,^ dv=H^'\x-x'). 

fca^co 7^^ J_i,_oo kaHi^^ {ka) 



V^i=+^ 2^£^^ m''{x-x')dx\ (57) 



The expression in (57) is the exact scaler function for the fields of a slot on a fiat sheet w^th a 
constant electric field VJ<:I2xq across a width of {2xo/k). For a sufficiently small Xq, (57) 
becomes 

^,i,^+^i?r(^) (58) 
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and for a sufficiently large x 



^to^ + 



iVoe si 



in Xq / 2 
•To V ^^ 



(59) 



Since (59) is the first term of the as3^mptotic expansion of (57) and (59) is also equal to 
(55), it is reasonable to assume that a more exact expression for (55) would be (57). Therefore, 
we can say that if Uq is the expression for a slot on a fiat conducting sheet, then the expression 
for a slot on a cylinder of large (ka) and values of (p small is 



*t=u. [(i-fv+n, »■+ . . :)-i (f »-4 ,'+i^, f+..)] 



(60) 



where 



y-- 



(kacp) 



3/2 



ka 



A plot of the bracketed term of (60) is shown in figure 3.7. This result corresponds exactly 
with the conclusions of J. R. Wait in his analysis of currents induced on large spheres [Wait, 
1956]. 

For a very narrow slot the expression of Uq should be, from (58), 



Uo~- 



■.+iY^m^\x) 



(61) 



where x is the distance from the slot. Our derivation for the corresponding term, shown in (55), 
has only the first term of the asymptotic expansion of the Hankel function appearing in (61). 
Therefore, for narrow slots (61) may be used in conjunction with (60). 

In conclusion we have, for large radius cylinders {ka> 10), the expression for the magnetic 
field on the surface of the cvlinder 



oofxe 



(62) 



where ypt aiid xp]^ are the traveling wave components in the positive and negative cp direction. 
The expressions for xj/ty ^or a slot located at <^=0°, has been discussed in detail in this section. 
From (35) we have 

^PA<p)=^l^n27r-<p) (63) 

which indicates the symmetric distribution of the magnetic field about the slot. 

Figures 3.8 and 3.9 show the plots of the ratio of the traveling wave component of the 
magnetic field due to a narrow slot on the cylinder to the field of a slot on a flat sheet. The 
computations are based on the approximations expressed in (35) and (54). 
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Figure .3.7. The curvature term amplitude and phase 
variation (after Wait [1956]). 
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Figure 3.8. Magnitude of the ratio of the traveling 
wave current on the cylinder to that of a flat surface. 
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Figure 3.9. Phase of the ratio of the traveling wave 
current on the cylinder to that of a flat surface. 

3.4. Comparison of the Harmonic Series and Residue Series Evaluations 

In the evaluation of the surface fields by the harmonic series, we were able to show the 
order of the errors involved in the final evaluation. However, in the residue series evaluation, 
we did not attempt to indicate the order of the errors involved in the approximations. How- 
ever, it is quite obvious that as ka increases the residue series evaluation becomes increasingly 
accurate. On the other hand, as ka increases, the use of the harmonic series becomes con- 
siderably more difficult. Since there is a range of values of ka, where both methods can be 
employed, we devote this section for a comparison of the results of the two methods. We 
perform this comparison by numerical computations for C3dnider sizes with ka=10, 15, and 25. 

Figure 3.10 shows the results of three computations. The solid curve is a plot of the mag- 
nitude of the magnetic field computed by the harmonic series (see sec. 3.2). According to the 
estimate of the remainder (see eq (24)) the order of the error by this method is (1/M) for small 
(p and 



ikb<^h'(Sr] 



for large M(p. Therefore, the errors in the computations of the solid curves of figure 3.10 are 
of order less than 1 percent over all values of cp. The points, in the figure, designated by the 
symbols and A were computed by the residue series method. For the set of points o, the 
computations were based on the first order approximation of the roots, Vm2, of the Hankel 
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Figure 3.10. Comparison of residue series and 
harmonic series evaluation. 
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functions discussed in appendix II. On the other hand, the set of points, A, were computed by 
using a second order approximation for the roots Vm2' It is quite apparent from the figure that 
the second set of points are closer to the harmonic series evaluation, as expected. Actual com- 
parison of the numerical values shows that the latter set is off by less than 2 percent of the har- 
monic series evaluation. A comparison of the phase of the field by the two methods, which is 
not shown, has a maximum deviation of 5 degrees for either set of roots. 

It should be noted that, in the figure, the region very close to the slot was omitted for con- 
venience. The agreement between the two methods of evaluation holds for points up to several 
slot widths. However, if the field distribution across the slot and its edge is desired, the har- 
monic series evaluation will yield a more accurate result for values of ka<C2b. The expression 
in (60) can be used for ka^lO provided the expression for a slot on a flat infinite sheet is known. 

4. Circumferential Slot Excitation 

In this section we will consider the case of a uniformly excited circumferential slot running 
entirely around the cylinder. 

Describing the field distribution on the surface of the cylinder we have, with F representing 
the voltage across the slot, 



From (2-10) we have 



Bn2 (f) =^ £' dZ jj g, {<p) g,{Z) e^^^^-i'^d<p. 
Using (1) in the expression above, we have 



(1) 



2tJ J_. rnr [^0,715^0. 



(2) 
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Substitution of (2) into (2-13) and (2-14) loads to 

and 

^^' 2:r }c,{k'-t')H^^KHk^-i') ^ 

where Ci is the contour defined in figure 2.2. 

Using (3), the equations in (2-8), and solving for the magnetic field on the surface of the 
cylinder, we have 

Hz{r=a)=0 
and 



Tlic ])i()l)lom, therefore, reduces to the evaluation of llic iiilciir-i 



'^'^=^iL, 



aVF-rH«'(aVF-f^) 



di. (5) 



For Z>0, we can deform the contour Ci about the branch cut as shown in figure 4.1. The 
contours Ci and C2 are equivalent since for i?i— >oo the integrand converges to zero and the 
integral about the outer circle does not contribute to the value of /. We will next show that 
we can obtain an asymptotic expansion of (5) for kZ<^(2kay. 

The argument of the Hankel functions in (5) varies between 7r/4 to — tt along the contour 
C2. Therefore, with a sufficiently large Bq, we can take the asymptotic expansion of the ratio 
of the Hankel functions appearing in (5). With y=a-JW^^^, we have 

m^'iy) \ '^iy 2{2iyy^{2ivf 8(2«/)^^"^'^ 'J 

Substitution of (6) in the integrand of (5) leads to 



(6) 



Therefore, in general, the asymptotic expansion of I{kZ) will have the form 

m—l 

I{kZ)=-^ a^g^{kZ) + ibJ„{kZ)+R^ (8) 

R^=0^{ka)-''^j^ JjUfyi^dij (9) 

and the first few coefficients a„ and b„ are: 

n 1 -1 



2(/:a)^ ' 8{kay 
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Figure 4.1. Equivalent contours for the integral I(Z), 






b,=- 



-25 1 



' 128 (kay 
From complex variable theory (10) becomes 

Residue L(wrl=- J='(^^=l)! ^T^ 1(1^^74=-.; 
From Courant and Hilbert [1953, page 482], we have 



^"(^^=-^^ "tw (f T I ^"^^(f^-i)^-''^^f- 



Therefore, (11) becomes 



t(-) 



(^^D 



Using the well known identity 



we finally have for ^11) 



H'X{kZ)=e-'^^m\kZ) 



(0 






(12) 
(13) 
(14) 

(15) 

(16) 
(17) 



It now remains to determine (9), which is the order of magnitude of the remainder term 
in (8). ForO<A:Z<l, 



and for fcZ> 1 



Therefore 



J 0,(1 -n™ L \2^vJ 



(18) 



{kaY 



Therefore (5) can now be expressed in the following asymptotic form with sufficiently 

large : 



354 



1 T- -ikZ 



e-i^z(l + ikZ) , 1 



mikay 



ka 



m^'ikZ) 



+1 (i-3 ^^'^^'^^^-ik w ^^^^ ^'^)]+^^ [ w^- ^^^^ 



Therefore the magnetic field on the surface of the cyHnder excited by a circiimferontial 
slot is, for ^«1 

rr (y. ^ m^ui^y. fi i { ^ l + ikZ \ e'^^' , 1 kZ Hj'^kZ) 25 (JcZl HlHJcZf] 

ii^[z.)= 2 ^0 ^'^^^L V^ita 2{2kay)H^'KkZy2{2kayW^\kZ) 24. {2kay W^\kZ) } 

(20) 

It is quite obvious that the bracketed term in (20) approaches unity as {ka) approaches 
infinity. As we let {ka) approach infinity, the problem reduces to the case of an infinite narrow 
slot, which is also the limiting case of (20). 

Figures 4.2 and 4.3 are plots of the magnitude and phase of 

^(7^-^^±.(^_l±^\ ^''"^ . ^ ^^ Ii'^'{kZ) 2^, {kZy m'^kZ) 
^^^^-^~^2kaV 2{;2kay)m'\kZ)^2 {2kay H^'^kZ) 24.{2ka)' m'\kZ) ^^'^ 

where various values of {ka) wei'e used as parameters and plotted as a function of normalized 
axial distances, (Z/2a). It can be seen from the curves that there is an increase in magnitude 
of /(Z) which shows the contribution of the various sections of the source whose path of prop- 
agation is along a spiral path on the surface. 
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Figure 4.3. The phase of the ratio of the magnetic- 
field on the cylinder to Ho^2)(k2). 



5. Mutual Admittance of Slots Coupled Externally 



In this analysis we wish to express the Mutual Admittance of a pair of slots coupled to each 
other by their External Fields. By ' 'External Fields^' we mean the electromagnetic field region 
excluding the transmission line regions feeding the slots. To obtain the expression for the 
^'External Mutual Admittance'^ we shall assume that the pair of slots are fed by independent 
transmission lines and that the slots are shunt elements to their respective feed lines. 

The external magnetic field of a single slot can be expressed as (see Oliner [1957]) 



HM=^^^(j^>^E)^rndS 



(1) 
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where E is the electric field across the slot, j^^ is the dyadic Green's function for the exterior 
region of the slot and S is the aperture of the slot cut in a conducting surface. Therefore, the 
external field due to a pair of slots can be expressed as 

He{r)={{ (^XE,).^ndS+{{ inXE2)rndS (2) 

where Ei and E2 represent the respective electric fields across the slots, Si and S2 their respective 
apertures. 

The fields of the slots in their transmission lines can be expressed in terms of voltages and 
currents. According to Marcuvitz and Sch winger [1951], the internal magnetic field due to a 
slot in a waveguide is expressed as 

Hi(r)=^ {I'''+r'')P^'\r)-i ^ (y^^>+F^2))^(2)(^)_^ JJ (^y^^) • ^dS (3) 

where /? is the Green's Function Diadic for the waveguide and ^^^\ /3^^^ are the standing wave 
type mode functions for the waveguide. 

The continuity of the tangential field required that Hi(r)=He(r) at the slot. Therefore, 
using (2) and (3) for slot (1) we have 

J£ (nX^i) • {n+i0)dS+jj^ ilxE,) .ndS=\ (/i"+/P)ft<»(r)-i|-"(Fr+Ff>)/3f' (r). 

(4) 

If this slot and its external load can be represented as a shunt element in the transmission line 
then 

/,:=7a)=_/(2) andFi=Fi(i>=Fi^2) 

thus reducing eq (4) into 

-iY,Vifi['KT)={{ inXE,)^{^n+imS+{{ Cnxh-rndS, (5) 

J J Si J J S2 

But, by the definitions given by Marcuvitz and Schwinger [1951] 

2lt=iYojj^ (nXEi) . ^{2) dS=I[ (6) 

where /^ is the total current in the shunt admittance. Multiplying (5) by (nXE) and integrat- 
ing over S, we obtain 

"^'^'^Sfs [/£ (nxA) ' {rn+iPi)ds'j^ {nXE,)dS 

+jj^\^jj^ {nXE2)'rnds'^-(nXEi)dS. (7) 

We can write a similar expression for the second slot in terms of its characteristic equations, thus 

_F,/;= (( r f( (nXE2) . {rn+i^2)ds]' {nXE2)dS 
J J S2L.J J S2 J 

"iX [JX (^X'^i)-^*'^^]- CnXE^)dS. (8) 



I Si 

+ 
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Figure 5.1 represents an equivalent circuit whose elements we shall identify in terms of the 
expressions discussed above. Writing the network equations for the equivalent circuit we find 






(9) 



Multiplying the above equations by F^ and V2S respectively and considering ideal transformers 
(1) and (2) we have 

VJ[=VuIu=Y'nVl+YnV2yu (10a) 



VJ'2=V2j2s=YnV^sV2,+Y',,V\ 



(10b) 



Comparing equations (10a) and (10b) with (7) and (8) we immediately observe the equiva- 
lence, which leads to 



Yr 



Vr.V^s 



(11) 



where Vu and F25 can be the slot voltages and the bracketed integral represents the magnetic 
field due to slot (1) across slot (2). 

Considering the short circuit condition of the secondary (1^2=0) and eqs (6), (7) and (10a) 
we obtain an expression for the self impedance of slot (1) at the transmission line terminals 



A^fJJsiLJJ^i 



^riXE{)-i;/n+iPi)dS [(nXE,)dS 



]■ 




(12) 



FiGUEE 5.1. Equivalent circuit of coupled slots. 

357 



Similarly the self impedance of the second slot is 

1 N' ^^S^^,^^^'>^'^''^^^'^^^'\^^^ 



■^ 22 



TiX ^^x^2)-^2''^^T 



(13) 



These are the same expressions as those obtained by Oliner in his study of impedance of 
narrow slots [Oliner, 1957]. It remains to show the expressions for Ni and N^. Considering 
the short circuit conditions again, we have 



where, 



Therefore, 



Vx-N^u or iVi= 



Similarly 



N:= V V ^14) 

-'^ 11^ Is 



N.= -^-^^ y (15) 

Physically, A^i and A^2 represent the conversion factors of the slot voltages to the trans- 
mission line voltages of the waveguides. Such a number is necessary to evaluate the coupling 
coefficient, which can only be measured from one of the transmission lines feeding the coupled 
system. 

sin n\pQ cos ncp 
,=1 n^po n^ 



to 

6. Appendix 1. Evaluation of 8^=^ 



In this appendix we evaluate the sum Sk of the infinite series 

^ " sin nxf/Q cos ncp ,^. 

for k=l, 2, 4, ... . With k odd, the series is rather difficult to evaluate. Since we have 
avoided these cases in the text, with the exception of k=l, we shall discuss the series with k 
even and then evaluate Si for values of \//q small. 

We can write (1) as the sum of two series as follows: 



S. 



_ J_ r^ sinnjip+xPo) ^ sinX^-^n ^ . 



We have from Jeffreys and Jeffreys [1956] for integer values of r^l 

hi n^'-' ~ 2(2r-l)! ^'^"^^^^ ^^^ 

where \p2r-i(t) are the Bernoulli polynomials and can be represented by the coefficients of 
Z^/r\ in the expansion 

-^^=^m^' (4) 
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Substituting the expression in (3) into (2), we obtain for even k>2 



S,=^ 



1 (27r)*+'(-l) 



f[,.„(^y,.„(^)} 



"= 2^„ 2(k+l) 
We evaluate the specific case of k = 2 and k=4. From (4) we have 

Ut)=t'-\t'-\-lt'-\t. (6) 

Substituting (6) into (5) and simplifying we finally obtain 

'^^=90-48+12 ^ -(,12+24; " +12 "''«^-V36+240r»- ^^^ 

We now proceed to evaluate Si, for which one cannot employ the identity in (8). We have 
g ^J_ y^ sin nypQ cos rup ^ 1 T^o T " cos njip+xf/o) ^ cos n(<;g— ^o) "j , 

^0 n = l ^ ^l/^o J Lw=l ^^ n=l ''^ J 

Now consider tlie following function for 0<arg a^<7r 

Then, for x real and — 7r<.r <7r, we have 

|:^^=Ee[5(x)] = -ln2(sin|). (10) 

Substituting (10) into (8), we obtain 

For i/'o small and \pi)<<p<iv — \pQ, we have 

^'-^J'" t'^ ^ (^"' l+t ^°^ 1)+^" ' (^^'^ i-t ^-^^ i)] ^^°- ^^') 

Carrying out the integration and simplifying we obtain 

o 1 ^^ aV ' 2^ /'AoV 2nitan<^/2, ["tan (<^/2)-(i/^o/2)"l . 1 ./x /-, ^x 

^,=.1 -.^- In 4 [sm^ --(--) cos^ 2]+-^ ^" [ tan (W2) + W2) J+2^o^-^^^ ^'^^ 

where A ((i?) is the constant of integration. From (8) and (10) we have 

Si(^,\^o=O) = -ln(2sin0 (14) 

From (13) we have 

lim ,5,(^,^0) = -In 2 sin |+lim-^- (15) 

Tliis will only be satisfied if /i(<^)=0. In (13) the third member of the expression at first 
glance seems to be singular for ^o=0. However, the limit of this term as ypQ-^Q is unity. 
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Therefore, (13) becomes 

S.^l+(^-r) ,„ (2 ,i„ l-fc CO, l)-(^+'^ In (2 s„, |+*. cos I). (16) 

Special care is needed to evaluate Si at the point where i/o=(p and <p=7r. When we take the 
limit of Si as (p approaches xpo and t we obtain 

Sl(^=^Po)=l-^n2^Po 

Sii<p=T) = -ln2. (17) 

To evaluate Si for 0<(p<\l/Q,we need to reconsider the problem beginning with (8). Thus 

Si=^± r [± ^osn(^+^o) ^^ cosn(^-^o) l ^ (,3) 

Using (10) to evaluate the sums in (18) we take for cp small and 0<(p<yl/Q 

Carrying out the integration and simplifying, we obtain 

*-^ [2 ..„ t >„ 4 (..„^ f -(|)>4 ..„ f +. 1„ gii||tg]. « 

As a check we note that at (p=\l/o (20) agrees with (17) if we assume ypo small as in (17). An 
important point to obtain from (20) is Si at (p=0. The last term in (20) approaches zero as 
<p approaches zero. Therefore, 

Si(V^o,^-0)=| tan I j^l-ln (2 tan |)]. (21) 

If ^0 is small (21) becomes 

>^i('Ao,^=0) = l-ln^o. (22) 

This completes the evaluation of Si for \f/Q small for all values of (p. 

7. Appendix II. An Expansion of Hankel Functions in Terms of Airy Integrals 

In this appendix we will show an expansion of the Hankel Function of large complex 
argument in terms of the Airy integrals and its derivatives. This expansion is generally 
employed for real values of the argument of the Hankel function [Bremmer, 1949; Wait, 1959]. 
However, it is easy to show that the expansion is valid for the half-plane defined by 

-(^+*)<argZ<^-. (1) 

where Z is the argument of the Hankel function. 

We consider the integral representation of the Hankel function 

^(2)(Z)=-i. r g-^Binhr+rr^f (2) 

'T^ Jli 

where the path of integration is as shown in figure II. 1. To investigate the region of con- 
vergence we take 

Z=x+ iy and f = J + iV' 
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Thou 

Z sinli ^=x sinh J cos r)—y cosh ^ sin 7]-{~i{x cosh f sm rj+y siiih ^ cos 77). 

The integral in (2) will converge if, along the patli of integration, 

Re [Z sinh ^] = x sinh ^ cos rj—y cosh { sin ?7>0. 

This requirement is satisfied by 

-^<arg Z+€<C^ 

where 

Therefore, we can choose an appropriate path L2 to define H^^^ (Z) for 



(^-^+cj<argZ<^-e 



(3) 



which defines the lialf plane of convergence (see fig. II. 1). 

We now use tlie Fock expansion [Fock, 1945] with the substitution 



/7\i/3 /7\i/3 

.=Z+(f) ^andX=(f) f 



and the assumption that t and X are finite and \Z\ is large. 
This expansion gives 



dx 



where F is an equivalent path of L2 in the X-plane as shown in figure II. 2. 
We now need to investigate the properties of the integral 



-Jir Jr 



X3 



known as the Airy Integral. 

The function w{t) satisfies the differential equation 

'w"it) = tw(t) 



(4) 



(5) 



(6) 



(7) 



Figure II.l. Path of integration and half-plane of 
convergence for the integral of the Hankel function. 
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Figure II. 2. Path of integration for the airy integral. 



X-PLANE 

which leads to the recursion formulas : 

w''(t) =2w'{t)+tw''(f)=2w'(t)+t^w{t) 
w'(t) =4:tw(t)+t^w'(t). 



We also have 



w 



it)-\Jl e-'i i-ty'H^i [I i-ty^'j 



(8) 
(9) 



where by {—t) we mean (te'""). 

Using the relationship in (6), (7), and (8), we have the asymptotic expansion 

The first term of the expansion in (10) represents the so-called '^Hankel Approximation^' 
and is used by several authors [Bremmer, 1949; Wait, 1956; Sensiper, 1957] in Uheir residue 
series evaluations. Using (9) and the first term of (10) we have, as Z-^oo 

mKZ)^^ dj" «K<)^(iy'y^ e-'"^H[% [I i.-tr^'] (11) 



where 



vl/3 

v-Z=(%) t. 



From (11) we have, for the derivative of the Hankel function 



But 



dZ= \z) 
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Therefore 

Hr{Z)^^(^"\\t). (12) 

Using (9) in (12) gives 

^-Z7r/6 /2\2/3 ro -] 

H«>'(Z)^^(^-|j </^«.[^|(_i)3/2j. (13) 

To complete this appendix, we investigate briefly the roots of H^^\Z)=0 and nf^\Z)=Q. 
These roots, to first order, can be obtained from the approximations discussed above. For Z 
large, snch as Z>50, the first order approximation for the roots is quite adequate. However, 
for moderate Z, higher order approximations might be necessary. Sensiper [1953 and 1957] 
discusses these roots quite thoroughly. At the end of this appendix we have, tabulated, the 
roots of H^^^' (Z)=0, to first and second order, for values of Z=10, 15, 20 and 25. 

From the approximations in (11) and (12), we can find these* roots, to a first order, quite 
easily. They correspond to the roots of 

where 

From Watson [1952] we have the identity 

Hl%m-zr;^, [J-vM-e'^JrMl (14) 

bill TT/O 

If we let ^=xe^'' and use the identity 

J,{xe''') = e''''J,{jr) (15) 

then (14) becomes 



te 



-IS 



^iH(^)=^^ [J-vzix)+Jvs(.x)]. (16) 

Therefore, to determine the roots Vna of Hf\Z) = ^ to the first order, we need to deter- 
mine the roots of 

J_l/3(x)+Ji/3W-0. (17) 

The roots Xmi of (17) are well known and we have a one to one correspondence between 

v,ay and Xr,ii. 

9 

[J — ^ ^iir — tL ( / ^3/2 

t^^=Q^y" e-'^i\ (18) 

The roots of H2/3(/3)=0 can be approximated in a similar manner. We have 

i2ir 



mm=-;\\j-2,M-e 3" j./3(^)]. 



(19) 
sm- 
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Letting fi—xe^'' again and using the identity in (15) we have for (19) 



Hi%{xe'^)-- 



%e 



2ir 

-i — 
3 



. 2t 

sm-r- 



[J-2/3{x)-J2/3{x)]. 



(20) 



Therefore, to determine the roots Vm2 of H^^^' (Z) to first order, we fiirst need to know the 
roots of 

J_2/3(^)-J2/3(^)=0. (21) 

The roots x,n2 of (21) are well known and we again have a one to one correspondence 
between Vm2 and Xm2- 



v„2=z+(^^y 



fim2 — Xjn2^^ 



-(I) i-i..y 



(22) 



O^ „2/3 

/ ^-^^2 -2 7r/3 

^m2 — r> *^ 

Sensiper [1957] in his discussion of these roots, gives the following expansions. 



(23) 
(24) 



Table II. 1 is a list of the roots, Xm and tm, and table II. 2 shows a comparison of the values 
of Vm2lZ computed by (22) and (24) for several values of Z. 

Table II. 1. Table of roots related to the roots of the Hankel functions 



m 


a: ml 


Xrn2 


tm:ei''3 


tn,2ei'fi 


1 


2.383 
5.512 
8.647 
11. 787 
14. 927 


0.686 
3.902 

7.058 
10. 202 
13. 348 


2.338 
4.089 

5. 358 
6.787 
7.944 


1.019 


2 


3.248 


3 


4.820 


4 


6. 163 


5 


7.372 







Table 11.2. Comparison of values of Vm2lZ computed by equations {22) and (24) 



m 


Vr„2lZfor Z=10 


Vm2/ZtromZ=15 




By eq (22) 


- 


Byeq(24) 


By eq (22) 


By eq (24) 


1 


1. 09-2*0. 161 
1. 28-20. 481 
1. 41-20. 714 
1. 53-20. 913 
1.63-21.09 


1 09-7*0 168 


1. 071-2*0. 123 
1.211-2*0.366 
1.313-2*0.543 
1. 401-2*0. 694 
1. 479-2*0. 830 


1. 069-2*0. 127 


2 


1. 27-?:0. 491 
1. 40-2*0. 693 
1. 51-20. 946 
1. 60-2*1. 14 


1. 208-2*0. 372 


3 


1.306-20.555 


4 


1 390-2*0 713 


5 


1. 466-20. 857 








Vm2lZ for Z=20 


Vm2lZ for Z=25 




Byeq(22) 


- 


By eq (24) 


By eq (22) 


By eq (24) 


1 _ 


1 059-20 102 


1. 051-2*0. 105 
1. 173-2*0. 308 
1.255-2*0.459 
1. 32.5-20. 589 
1.387-20.707 


1. 051-20. 0878 
1.151-2*0.262 
1. 224-2*0. 388 
1. 287-20. 496 
1. 343-2*0. 594 


1 050-?*0 0895 


2 


1. 175-iO. 304 
1. 260-2*0. 451 
1. 333-2*0. 576 
1. 398-20. 689 


1. 149-2*0. 265 


3 


1. 220-20. 394 


4 


1. 281-2*0. 506 


5 


1. 335-/0. 608 
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